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Abstract 

The role of reduced dimensionality and of the surface on electron-phonon (e-ph) coupling 
in silicon nanowires is determined from first principles. Surface termination and chemistry is 

o : 

found to have a relatively small influence, whereas reduced dimensionality fundamentally al- 
ters the behavior of deformation potentials. As a consequence, electron coupling to 'breathing 

> 

Q\ ' modes' emerges that can not be described by conventional treatments of e-ph coupling. The 

C ■ ■ consequences for physical properties such as scattering lengths and mobilities is significant: 

o 

fsj . the mobilities for [110] grown wires is 6 times larger than for [100] wires, an effect that can 

o . 

not be predicted without the form we find for Si nanowire deformation potentials. 



Silicon nanowires, beyond having been successfully demonstrated as conventional semicon- 
ductor devices," are beginning to be seen as the important blocks for novel energy harvesting 
applications, such as solar cells^ and efficient thermoelectric devices.- 1 ^ All these applications 
rely on a high electronic conductivity of electrons, while in the latter a low phonon conductiv- 
ity is also essential.-*"* At the heart of understanding the interrelation between these properties is 
the interaction of phonons and electrons: scattering of electrons with phonons reduces the con- 
ductivity of these devices, increasing heating and reducing their efficiency, while the scattering 
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of phonons with electrons reduces thermal conductivity, yielding better thermoelectric attributes. 
But although electron-phonon coupling in bulk semiconductors has been extensively studied, the 
nanoscale regime is not so well understood. 

Materials when structured as low-dimensional systems are known to exhibit new behavior rel- 
ative to their macroscopic properties. In the case of silicon nanowires, quasi-one dimensional 
structures, band folding can result in an electronic structure with a direct band gap. The Increased 
surface to volume ratio permits the surface chemistry to significantly alter the value of band gaps 
.— »^ In this paper, we demonstrate that reduced dimensionality in silicon nanowires also has a 
significant impact on electron-phonon (e-ph) interactions that has not been previously anticipated. 
We derive a new form for deformation potentials in Si nanowires that enables a general approach to 
the calculation of e-ph couplings, an approach that may be readily extended to other semiconduc- 
tor nanowires. The validity of such macroscopic theories at the nanoscale is central to advancing 
semiconductor technology (see, e.g., the developments in the effective mass approximation— - —). 
To our knowledge, this is the first paper to address the effects of dimensionality on the deformation 
potential theory. 

For this purpose, we employ density functional theory (DFT) to calculate the band structure 
of Si nanowires and apply strain to extract the deformation potentials. DFT is known to give 
good agreement with experiments of the band gap pressure dependence 16 ' 17 and to give accurate 
deformation potentials that reproduce the n-type carrier mobility in bulk. 18 Focusing on electrons, 
we find that contrary to common assumption, the tabulated deformation potentials not only vary as 
a function of size and orientation but also vary with respect to the direction of the applied strain. 
This has a direct influence on the strength of the scattering of electrons by phonons travelling in 
different directions. We develop a simple theory to take into account the full anisotropy of the 
deformation potentials and its impact on electron mobility. The effect of the direction of growth 
of the wires and surface passivation is also studied. Notably, the deformation potentials in [110]- 
oriented wires are found to be highly anisotropic when compared to those of [100] wires or bulk 
silicon. This results in the suppression of the scattering from breathing modes, and, coupled to 
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their lower effective mass, leads to much higher mobilities for [1 10] wires. 

Electron-phonon coupling in nanowires can be expressed in a similar manner to the deforma- 
tion potential theory for anisotropic scattering in multi-valley semiconductors.— A first approx- 
imation might be to use the bulk values for the deformation potentials to determine an isotropic 
model independent of the nanowire diameter or surface termination together with the effective 
masses calculated for specific nanowires. Unfortunately, as Leu et al— find in their recent paper, 
the deformation potentials of the valence and conduction bands of [1 10] and [1 1 1] Si nanowires are 
substantially different from those in bulk Si for strains applied along the nanowire axis. A common 
practice is to assume an isotropic e-ph scattering rate with an effective value to account for defor- 
mation potential variations in nanowires relative to bulk,— - — greatly simplifying the calculation 
of physical quantities derived from the deformation potential. Speaking against this approach is 
the fact that Fischetti and Laux 24 have shown that the anisotropy of the e-ph matrix element has a 
substantial effect on the magnitude of the electron-phonon scattering rate in a 2D electron gas in 
silicon. 

The effect of an acoustic phonon on the electronic band structure is equivalent to that produced 
by a slowly varying strain in the direction of the phonon. This strain introduces a correspondent 
slowly varying potential that scatters the electrons. The change in the energy of a band per unit of 
strain is called a deformation potential. The deformation potential theory for bulk semiconductors 
has been developed by Herring and Vogt,— who find the anisotropic electron- acoustic phonon 
matrix element for an elliptical valley in silicon to be: 

// q C (r) = {s rf 5R ? -q + S l( (k«-5R^)(q-k a )}cos(q-r), (1) 

where and Z u are the dilatation and the uniaxial- shear deformation potentials, respectively, 
q is the phonon momentum, k a is a unit vector parallel to the k-vector of valley a and oRq is 
the amplitude of the displacement of the atoms due to a phonon with polarization £. We will 
initially assume that this theory still holds for SiNWs, and find and E u from electronic structure 
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calculations. Our assumption will hold for wires with large enough diameter D. For quasi-one 
dimensional nanowires (D < 3 nm), 25 we show that corrections are needed. It will be seen that, 
for the quasi-one dimensional wires, becomes anisotropic with respect to the direction of the 
applied strain (or phonon); this effect is particularly pronounced for [1 10]-SiNWs. 

The band structure of nanowires can be understood to be due to the spatial confinement of 
the bulk electronic bands. For a [100] nanowire, the indirect conduction band (CB) minimum at 
±0.30f is formed by the bulk Si valleys at Ih/ciq [±£ , 0, 0], with energy E x , where £ = 0.85. The 
CB minima at T are composed from the four A valleys at 2%/ao [0, ±£ , 0] and Ik/ciq [0, 0, ±£], with 
energies E y and E z , respectively. From explicit band calculations for a nanowire, the four valleys 
that fold into T are not exactly degenerate owing to level splitting introduced by the surface. We 
can use a simple theory in which the effect of the surface is to introduce interactions between four 
degenerate bands, so that the CB minima in the SiNW can be described by the eigenvalues of the 
following matrix, 

1 E y a p 
a E y B 

E = 

P P Ez 
V P P a 

where a and /3 are the g-type and f-type matrix elements between the bulk- like bands, respectively, 
introduced by the presence of the surface. The four eigenvalues of this matrix are 
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E z + E y + 2a ± a/ 16/3 2 + (E z - E 



(3) 
(4) 



In the absence of strain, E x , E y and E z are degenerate and E\ = E2, which we find to hold for 
hydrogenated SiNWs. 26 Similarly, in a [1 10] nanowire, the conduction band minimum is at T and 
formed from the bulk valleys at 2n/ao [0, 0, ±£j], with energy E z and the band minimum close to 
±0.9^, is formed by the bulk valleys 2tz/qq [±£,0,0] and 2n /clq [0, ±£,0], with energies E x and 
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Ey. The effective mass of the lowest conduction band in the wires is thus given by the transverse 
mass in bulk, while the subband separation depends on the longitudinal mass. The effective masses 
are directly calculated for the wires from first principles (see | [table] [1][]1| ). The meaning of the 
effective masses m t and m/ for the wires is as follows: m t is the effective mass along the direction of 
the wire, since, from an effective mass perspective, the lowest conduction band can be considered 
as formed by the valleys with transverse mass along the direction of the wire. Likewise, m/ is such 
that it gives the correct separation between the first sub-bands. 



The valleys shift according to Eq. ( |[equation][l][]l[ ) as strain is applied, yielding the deforma- 
tion potentials for each of the NW bands. From these shifts and Eqns. © and ©, we can obtain 
Erf and E u . For example, a strain in the [100] direction yields, E^ = (E\ + E2 + E2 + E^/4 where 
Ej are the deformation potentials corresponding to the £, bands. The variation of a and /3 due 
to strain is found to be small and is ignored in the following. The resulting values for the wires 
studied in this work are she 
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Figure 1: Diameter and surface passivation dependence of the room temperature mobility of 
SiNWs as indicated. The dash and dot-dash lines show the the effect of disregarding the scattering 
by breathing modes in [110]- and [100] -oriented wires, respectively. 



We have calculated 27 the deformation potentials for the first conduction bands for wires of 
different sizes with orientations [100] and [110], by varying the interatomic distance by 0.1% to 



0.5% along and perpendicular to the wire axis. |[table][l][]l| shows the E^ and E u calculated using 
the values from the lowest bands at Y. Here, E^m and S^_|_ are the equivalent of the E c [ deformation 
potential in the parallel and perpendicular directions to the wire axis, respectively. From these 
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Table 1 : Calculated effective masses and splitting matrix elements for the conduction bands of the 
different SiNWs. Also shown are the bulk-like deformation potentials for the first CB of SiNWs. 
Zj± and S^ii are the Zj deformation potentials for phonons in the parallel and perpendicular direc- 
tion to the wire, respectively (see text). The brackets next to the atoms per simulation cell indicate 
silicon crystal orientation along the nanowire axis. For the [1 10] wires, a gives the splitting of the 
two lowest conduction bands. m e is the electron mass. 



Width [nm] 


Atoms/cell 


m t 




mi[m e ) 


a[eV] 


j8[eV] 


Srfll [eV] 


Z d± [eV] 


S ;( [eV] 


0.82 


S* 21 # 20 [100] 


0.46 


1.4 


-0.3 


0.038 


-4.2 


-4.2 


8.2 


0.77 


S/25#20[100] 


0.37 


0.9 


-0.039 


0.009 


-6.1 


-6.1 


8.7 


1.09 


&37#28[100] 


0.44 


1.1 


0.078 


0.025 


-4.7 


-4.7 


8.6 


1.26 


Si 51 H 36 [100] 


0.4 


1.1 


-0.074 


0.017 


-5.5 


-5.5 


9.2 


1.9 


S*9 7 // 4 4[100] 


0.39 


1.1 


-0.047 


0.0026 


-6.6 


-6.6 


9.2 


1.086 


Si 31 (OH) 2S [100] 


0.72 


0.95 


0.24 


0.05 


-4.8 


-6.5 


5.3 


0.53 


S/ 8 // 12 [110] 


0.15 


1.5 


0.58 




-9.4 


-1.3 


3.7 


1.09 


si 23 H 20 [no] 


0.15 


0.62 


0.38 




-10.0 


-1.8 


3.7 


1.15 


si 24 H l6 [no] 


0.15 


0.62 


0.39 




-9.1 


-0.1 


3.9 


1.35 


S/42#20[H0] 


0.15 


0.62 


0.21 




-10.5 


-1.3 


4.1 


3.07 


Si l52 H 40 [llO] 


0.15 


0.62 


0.018 




-9.5 


-3.3 


5.8 


1.15 


Si 24 (OH) l6 [ll0] 


0.19 


0.95 


0.03 




-7.5 


-3.2 


3.1 


1.35 


S142(0#)2O[11O] 


0.15 


0.95 


0.33 




-7.3 


-5.1 


3.9 



direct band structure calculations, we have found that Eq. ( [equation] [1][]1 ) still holds for small 
nanowires, except for hydroxylated and [110] wires, where the expression can be corrected to 
include the anisotropy due to the breaking of bulk symmetry, as seen below. Surprisingly, there is 
very little variation in the Zj and S M bulk- like deformation potentials with the strain direction in the 
[100] wires, and Z u has the same value as bulk silicon. 18 On the other hand, Zj varies substantially 
with strain direction in the [110] wires, and Z u is half its value in bulk. The anisotropy in Zj can be 
accounted for in Eq. ( [equation] [!][]![ ), by interpolating this deformation potential as a function of 
strain direction as follows, 



TtHO/^ _ n 2^UQ , / _2 , 2\ -.110 

W - Qx^d\\ + {1y + 1z ) ^di 



(5) 



where the qi are the cartesian components of the dimensionless unit vector q along the phonon 
momentum. 



In what follows, we use the values of | [table] [!][]! | and the Boltzmann transport equation to 
estimate the mobility of the SiNWs. The momentum relaxation time, Tk, due to acoustic phonon 
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scattering can be expressed as 



4 = > f d M-^) 
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X 



(6) 



where the upper and lower symbols in "±" indicate the emission or absorption of a phonon, L, 
are the dimensions of a wire of rectangular cross section, with x the periodic direction, is the 
velocity of carriers with momentum k, eq the polarization of a phonon with momentum q in branch 
£, and 



Gm Zl m y ,n' z ,n' 



X 



L z L y JQ 



Ly 



o 



dzdys'm(k z z) s'm(k y y) x 



sm(— n z z)sin(— n y y)e 

L Z Ly 



(7) 



25 



for the meaning of the indices and vectors 



is the overlap factor of the wavefunctions (see Note 
in this context). The choice of "effective mass" wavefunctions for Eq. © over the Kohn-Sham 
wavefunctions is based on ease of use and to simplify the understanding of the different contribu- 
tions to the scattering from the different modes. This approximation should work for large wires. 
In smaller wires, the electron confinement and the matrix elements introduced by the surface push 
five of the six sub-bands up in energy, beyond the range of acoustic phonon scattering at room 
temperature. As a result, this approximation works well for the lowest sub-band, which is the main 
contributor the the mobility. Some corrections should be expected for the scattering into higher 
sub-bands, but they are not important for the cases considered in this work. We therefore expect 
the main difference in comparing the different wires to come from the deformation potentials. Eq. 
© can be solved numerically, and using the elastic approximation does not yield noticeable differ- 
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ences. The acoustic phonon-limited mobility can then be obtained as fi(Ep) = „f i where m t 
is the effective mass of a carrier in the conduction band. Thanks to the breaking of the degeneracy 
of the conduction bands due to the surface matrix elements and confinement observed, we expect 
the contribution from optical and intervalley to be small for the range of energies considered in this 
work. We would like to point out that, while we have used the bulk phonon dispersions, discretised 
perpendicular to the wire axis, the effect of the full deformation potential remains the same, and 
we do not expect a large difference in the mobilities between this approach and using the phonon 
band structure of the wires.— 

The room temperature mobilities yl(Ep = 0) of wires with varying orientation and surface 



termination versus wire size are summarised in [figure] [1][]1 where it is seen that the nanowire 
orientation has a much larger effect on electron mobility than surface passivation. Due to the 
increased effective mass, OH surface passivation affects only the [100] growth orientation, halving 
the mobility compared to H-passivation for the case considered. A lower mobility is observed for 
the [100] -oriented wires compared to the [110]-SiNWs. One factor yielding the lower mobilities 



for [100] wires is their larger electron effective masses (see |[table] [!][]![ )■ However, the electron 
relaxation times in the [100] and [110] wires differ fundamentally in that the e-ph coupling for 
breathing vibrational modes in the [110] nanowires are suppressed relative to the [100] wires or 
bulk silicon. For the [1 10] wires, electron scattering due to the phonon modes along the wire axis 
is usually more than four times as strong as that due to the modes perpendicular to the wire axis 
(breathing modes), and larger than for the [100] wires. In the latter, however, the scattering is 
more isotropic, with the result that breathing modes have similar e-ph coupling strengths as for 
excitation of modes along the wire axis. To demonstrate the importance of the breathing modes in 



suppressing the [100] -mobility, in [figure] [1][]1 we consider the effect of de-coupling the breathing 
modes from the scattering matrix. As the nanowire diameter is increased, coupling to the breathing 
modes play an increasingly important role in determining the electron mobility. 

The anisotropy in the deformation potentials is also expected to influence the thermoelectric 
performance of NWs. The electrical conductivity and the electronic and phonon thermal conduc- 
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tivity that determine the figure of merit (ZT = ^ T, where Kj and K e are the ionic and electronic 
thermal conductivities, respectively, o is the electronic conductivity, S the Seebeck coefficient, and 
T is the temperature), will depend on several factors such as temperature, impurity concentration, 
surface roughness and chemistry. Considering the effect of the anisotropic potentials alone, better 
performance is expected for the [110] NWs compared to the [100] or bulk silicon: the enhanced 
scattering of the phonons by electrons along the direction of the wires (i.e. the direction of the 
thermal gradient) decreases thermal conductivity, while the overall reduction of the scattering of 
electrons by breathing modes enhances the electronic mobility. Hence, these two effects combined 
should result in a higher figure of merit for the [1 10] NW. 

We would like highlight the consequences of attempting to apply bulk-like deformation po- 



tential models to calculate the electron mobility in SiNWs. In [figure] [2] [] 2, we show the room 
temperature mobility versus the Fermi energy for a wire of 3 nm diameter oriented along a Si 
[110] direction, computed using the deformation potentials calculated in this work and compared 
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for bulk silicon. 



to the mobility calculated using = —0.029 eV and Z u — 8.77 eV, from Ref. 
The mobility obtained by using bulk deformation potential values underestimates the scattering 
by longitudinal phonons along the wire direction resulting in a significant overestimation of the 
electron mobility. It has been argued 22 that it is possible to build an approximation by which the 
value of an isotropic deformation potential reproduces the universal mobility curve obtained from 
wires of large diameter. However, this approach underestimates the electron mobility relative to the 
use of the explicit model developed in this work, as displayed by the green curve (labeled Iso) in 



[figure] [2] [] 2 Alternatively, to determine whether we can apply an isotropic model to reproduce 



the calculated mobility, we may attempt to define an effective isotropic deformation potential 



^effq 



(SR^q)cos(q-r) (8) 



that, used with an isotropic version of Eq. © (with the longitudinal phonon branch only) repro- 
duces the anisotropic mobility at Ep <C 0. The mobility obtained in this manner is also plotted 
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in |[figure][2][]2[ We find that the isotropic model is still not adequate in this case, where Z d is 



highly anisotropic, and Z u is smaller or near S^. Nevertheless, an isotropic model will work when 
the condition ~ |S M is satisfied, as is the case of the silicon nanowires oriented along a [100] 
direction. 



0.8 



> 

c l 

E 
U 

cn 
O 

a. 



0.6 



0.4 



0.2 





/ \ / \ 1 








— Full 

— Iso 








— - Bulk 








— Iso eff 






7 ,~\T \ \ 








1 \ K_ 








\ 1 

' ^—x \ 








V 



























0.2 



0.4 0.6 0.! 
E F [eV] 



1.2 



Figure 2: (colour online) Comparison of the room temperature mobilities vs Fermi energy of a 
3nm [1 10] wire between isotropic models (green, thin black), bulk deformation potentials (dashed, 
scaled down threefold to fit the figure) and the full anisotropic electron-phonon coupling (thick 
black). The bulk deformation potentials are those from Ref. 18, the isotropic deformation poten- 
tial is that used in Ref. [22I The thin black line represents the isotropic model using the effective 
deformation potential that fits the anisotropic model at Ef << (See Text). Interestingly, quan- 
tization effects corresponding to quasi-one dimensional electron density of states persists even at 
room temperature. 



In summary, we have calculated from first-principles the deformation potential of nanowires 
of different diameters, with [100] and [110] orientations and for hydrogen and hydroxyl surface 
termination, and extended the deformation potential theory to develop for the first time a model 
that incorporates the effects of lower dimensionality in the electron phonon scattering rate. While 
we find that surface termination changes the bands substantially, 1 ^*^ it has little effect on the de- 
formation potentials. To the contrary, we find that the deformation potentials in [110] are highly 
anisotropic, in contrast to those of [100] wires and bulk silicon. This effect results in the suppres- 
sion of the scattering from breathing modes in [110] wires, which, added to their lower effective 
mass, leads to a greater n-type carrier mobility compared to [100] -oriented wires. 

This material is based upon works supported by the Science Foundation Ireland under Grant 
No. 06/IN. 1/1857, and the European Union project Nanosil. We would like to thank Stephen Fahy 
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